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Abstract 

, The macroscopic dielectric function in the random-phase-approximation 

o . 

<^> without local field effect has been implemented using the local density approx 

in 

imation with an all electron, full-potential linear muffin-tin orbital basis-set. 
This method is used to investigate the optical properties of the semiconduc- 
tors Si, Ge, and GaAs under hydrostatic pressure. The pressure dependence 
of the effective dielectric function is compared to the experimental data of 
q | Gohi and coworkers, and an excellent agreement is found when the so called 



"scissors-operator" shift (SOS) is used to account for the correct band gap at 



r. The effect of the 3d semi-core states in the interband transitions hardly 
changes the static dielectric function, e M ; however, their contribution to the 
intensity of absorption for higher photon energies is substantial. The spin- 
orbit coupling has a significant effect on 6oo of Ge and GaAs, but not of Si. 
The Ei peak in the dynamical dielectric function is strongly underestimated 
for Si, but only slightly for Ge and GaAs, suggesting that excitonic effects 

might be important only for Si. 
71.10,+x, 71.25.Tn, 78.20.-e 
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I. INTRODUCTION 



The experimental determination of the optical properties of bulk semiconductors can 
now be obtained with high precision 0-0], yet our theoretical understanding is far from 
complete. The static dielectric constant, which can be obtained from a functional derivative 
of the electron density with respect to the total Kohn-Sham potential evaluated at the ground 
state, hence a ground state property, is over estimated by the local-density approximation 
(LDA) calculation ([§-[§]. The inclusion of the gradient correction to the pseudopotential LDA 
reduces slightly the discrepancy in the case of Silicon J7|. The underestimation of the E\ 
peak and the overestimation of the E 2 peak of the imaginary part of the dielectric function, 
e 2 (u), by one-electron band theory have generated theoretical work for almost two decades 
to account for these discrepancies. It was clear from the beginning that including excitonic 
effects, which have been detected experimentally |§ , could remove some of the disagreement 
with experiment p|-[T2|. However, the model calculations used to correct e 2 (u) have produced 
only a qualitative understanding of the problem. In particular, the latest model by Hanke, 
Mattausch and Strinati based on the time-dependent screened Hartree-Fock approximation 
and including both the local field and the excitonic effects, described correctly the E\ peak 
but underestimated significantly the E 2 peak of Si. The reason for the underestimation of 
E 2 was attributed to a bad representation of the band structure of Si by their Slater-Koster 
parameterization ||11|| . Del Castello-Mussot and Sham [12|] based their latest calculation on a 
k.p model around the L point (where E\ originates) and a multiple plane- wave model around 
the X points (where E 2 originates), and solved the Bethe-Salpeter equation containing the 
excitonic effect. Their model is an improvement over the non- interacting approximation: 
the E\ peak becomes stronger and the E 2 peak weaker. This model is very promising, but, 
being based on a k.p approximation to the band structure, it provided only a qualitative 
correction to the intensities of the E\ and E 2 structures. Calculations ignoring excitonic 
effects, but including the local-field effect, underestimated both E\ and E 2 peaks [[L3|]. 

One way to make theoretical progress in this field is to determine the correct contribution 
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of the one-electron theory to the optical properties of semiconductors. This allows us to 
define precisely the size of the many body corrections to the one electron theory. However, 
a common belief these days is that the eigenvalues and vectors of the Kohn-Sham (KS) 



equations have no direct physical meaning and hence should not be used to calculate 
optical spectra of materials. Only ground state properties derived from the total energy as 
a function of the electron density have in principle a direct physical meaning. 

While LDA was indeed intended to calculate ground state properties it could also be 
viewed as a simplified quasi-particle (QP) theory where the self-energy is local and static 
(S(r, r',t) w V xc (r)5(r — r')5(t), here V xc (r) is the local exchange and correlation potential 



as, for example, parameterized by Von Barth and Hedin p!5| . The KS eigenvalues are then 
QP energies and could be compared to experimental data. This argument is supported by 
calculations using the GW approximation of Hedin ||16||. These calculations showed that the 



valence QP energies of semiconductors are in good agreement with LDA and the conduction 



QP energies differ by approximately a rigid energy shift [|T7|,[T8|. In the literature this shift 
is often called "scissors-operator" shift (SOS) |J. 

First-principles local density approximation calculations started more than two decades 
ago, but the major problem of LDA, beside the well understood energy band gap problem 
19fl , is the numerical difficulty to determine selfconsistent electronic structure and optical 
matrix elements using a complete basis-set. The early ab-initio calculation of the optical 
properties of semiconductors by Wang and Klein, using a selfconsistent linear combination of 



gaussian orbitals produced static dielectric functions in good agreement with experiment |20 
But this agreement is fortuitous because the band gaps produced by this method are much 
larger than the LDA band gaps. The recent calculations of 18 semiconductors by Huang 
and Ching using an orthogonalized linear combination of atomic orbitals method produced 
LDA static dielectric functions which are in general smaller than experiment despite the 
fact that their band gaps are much larger than the all-electron or pseudopotential LDA 
band gaps pifl . Those underestimated static dielectric constants are most likely due to the 



incompleteness of the basis-set used in their calculations. 



Most of the theoretical studies of the optical properties of semiconductors in the literature 
use several approximations within the LDA, ranging from the use of spherical potentials || 
to the use of pseudopotentials || instead of all electron LDA potentials. In this paper we 
report precise calculations of the optical properties of bulk semiconductors Si, Ge, and GaAs 
under hydrostatic pressure using an all electron LDA linear muffin-tin orbital basis-set []22| . 



in which no shape approximation is made for either the potential or the charge density 
23fl . The semi-core 3d of Ga and Ge are included in a fully hybridizing valence basis set, 



and the rest of the core states are allowed to relax selfconsistently. The effect of spin-orbit 
coupling is also investigated. A systematic check of the /-sum rule is performed for all the 
calculations. We hope that this accurate LDA calculation will provide an excellent starting 
point for the determination of the local field and the excitonic effects in the optical spectra 
of semiconductors P-fl3|l. 



We have found that the static dielectric function, eoo, which is a ground state property, 
is overestimated by LDA over all pressure range, and that an excellent agreement with the 
experimental results of Goni and coworkers [TJ for 6^ of GaAs and Ge under hydrostatic 
pressure is achieved only when the so-called scissors-operator shift is used to account for 
the correct band gap at T. The inclusion of the 3d semi-core states of Ge and GaAs in 
the interband transition has almost no significant effects in e OT ; however, the 3d interband 
transitions contribute significantly to the magnitude of 62(0*) above 25 eV for Ge, and above 
12 eV for GaAs. The spin-orbit coupling increases the LDA values by about few percents. 

The rest of the paper is organized as follows. In section II, we describe the method of 
calculation of electronic structure and the macroscopic dielectric function based on our all- 
electron full-potential LMTO basis-set. In section III we present the electronic properties of 
Si, Ge, and GaAs and compare them to existing theoretical results. The calculated dielectric 
functions and a discussion about including the semi-core states and the spin-orbit coupling 
will be presented in section IV. In the same section we also compare our static dielectric 
function under hydrostatic pressure with the experimental results of Goni et al. [|TJ. The 
conclusion is given in Sec. V. 



II. METHOD OF CALCULATION 



A. All electron full-potential wave function 

The full-potential linear muffin-tin orbital method in its scalar-relativistic and full- 
relativistic forms |S| is used here to calculate the electronic structure and the optical prop- 



erties of Si, Ge, and GaAs under hydrostatic pressure. The Kohn-Sham [|1J] equations are 
solved for a general potential without any shape approximation p3| . In this subsection we 
describe the Bloch wave function inside the so-called muffin-tin spheres and the interstititial 
region. A correct determination of the crystal wave function is necessary for the accurate 
determination of the optical matrix elements. 

As for the cellular methods, the space is divided into non-overlapping muffin-tin spheres 
surrounding atomic sites where the Schrodinger or the Dirac equation for each principle 
quantum number v and momentum channel i is solved for a fixed energy E U £. In these 
muffin-tin spheres the trial wave function is linearized in terms of the solution of Schrodinger 
equation <p T £ and its energy derivative <p T e for the energy E U £, and for an atom of type r and 



momentum channel £ 22.24 



It can be shown that the Bloch wave function of site r calculated at site r' in the unit 



cell of the crystal at R = is given by [23 



Xr£m(r) |r' = fat'm' ~ T ) B^'^K , k) + ^fjm' - T ') B^^K , k) (1) 

I'm' 

Where B^, T im (K, k) and B^, T im (K, k) are renormalized structure constants obtained from 
the crystal structure constants By m , (m (K, k) to ensure that the Bloch wave function is 
continuous and differentiable at the boundary of each muffin-tin sphere. 

In the interstitial region, the muffin-tin orbitals are spherical wave solutions Hi to the 
Helmholtz equation with non-zero kinetic energy; these bases are Hankel functions for neg- 
ative kinetic energies or Neumann functions for positive kinetic energies k 2 . Such that each 
partial wave inside the muffin-tin sphere is allowed to have different kinetic energy, k 2 , in 
the interstitial region. In this region the Bloch wave function is given by: 



XreJr) = E e<kR ^(«» l r " T ~ RD^-(r — r — R) (2) 

R 

The interstitial-region Bloch function is expressed in plane waves over the reciprocal lattice 
using Fourier transform : 

X^ m (r) = E/K(k + G)e^ k+G ) r (3) 

G 

Where K = {r, £, m, Eg, k}, here the parameter E e is the linearization energy of the wave 
function in the muffin-tin sphere for the t momentum channel, m is the azumutal quantum 
number, and k is the variational parameter whose square is the kinetic energy in the inter- 
stitial region. The Fourier coefficients, /k, are obtained from a pseudo-wave function that 
is equal to the crystal wave function in the interstitial region and represented by a smooth 
function inside the muffin-tin spheres. The exact shape of these pseudo-functions inside the 
muffin-tin spheres is not important. The only requirement is that they are continuous and 
differentiable at the sphere boundary and have zero slop at the origin of each sphere. The 
plane wave expansion is multiplied by a three dimensional step function so that the wave 
function is kept only in the interstitial region. The knowledge of the Bloch wave function in 
the whole unit cell allows us to calculate the Hamiltonian and overlap matrix elements in 
order to solve the effective one electron Schrodinger equation. 

Three different kinetic energies were used for each subset of s and p derived bases in the 
basis set; two kinetic energies were used for bases derived from orbital parameters £ > 1. The 
basis sets used in calculating total energies and structural properties were for Si: 3(3s3p), 
2(3(f), for Ge: 2(3d), 3(4s4p), 2(4d), and for GaAs: 2(Ga 3d), 3(Ga AsAp), 2(Ga Ad), 3(As 
AsAp), 2(As 4c?) ; the pre-multiplicities in this notation refer to the number kinetic energies 
used in this basis subset. The basis functions for each material comprised a single, fully 
hybridizing basis set. Note the presence of both 3c? and 4c? derived bases on Ga and Ge. 
A useful feature of the method used in these calculations is the ability to incorporate basis 
functions derived from the same orbital atomic quantum numbers but different principal 
atomic quantum numbers in a single fully hybridizing basis set. This feature entails the use 
of multiple sets of radial functions to represent bases with different principle atomic quantum 

6 



numbers. This capability was particularly useful in calculating the high-lying energy bands 
which were used to obtain the dielectric functions to high energy; the basis sets employed for 
this purpose are given in Table I. Seven to eight kinetic energies were used in the basis sets. 
Accurate resolution of the bands to high energy was necessary to converge the calculation of 
the real part of the dielectric function, which was obtained from the imaginary part through 
the Kramers- Kronig relation. An interesting consequence of the relaxation of the Ga 3d 
states as valence states is a significant decrease in the calculated band gap [p5 |. 

For the core charge density, the Dirac equation is solved selfconsistently, e.g., no frozen 
core approximation is used. The exchange and correlation potential is treated within the 
Von Barth and Hedin parameterization [|1J. To account for the relativistic effects in the 
dielectric function, the full-selfconsistent relativistic band structure is produced by including 
the spin-orbit coupling to the Hamiltonian. In Table I we show the orbitals used to describe 
the electronic states of Si, Ge, and GaAs. This large number of orbitals is necessary to 
calculate accurately the eigenvalues and eigenvectors up to 5 Ry above the highest valence 
states. These electronic states will be needed to determine the dynamical dielectric function 
and the converged static dielectric function through the use of Kramers- Kronig relations. 

The completeness of basis-set, with different variational k values for each partial wave in 
the interstitial region together with the Fourier representation allows the method to treat 
open structures like the zinc-blende structure studied here without having to resort to the 
so-called empty spheres | 26|| . The high energy states are also determined more accurately 
due to the use of many k values. As a test we show in Table II the eigenvalues of Si at high 
symmetry points of the Brillouin zone compared with some recent results from first principle 
calculations based on ab-initio pseudopotential and Gaussian orbital methods p7| , |28[| . The 
agreement of our calculation with the previous calculations is excellent. 
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B. Dielectric Function 



Here we give a concise review of the determination of the dielectric function of a semi- 
conductor crystal due to the application of an electric field. We also determine the approx- 
imations used to obtain numerical results for Si, Ge, and GaAs under hydrostatic pressure 
with or without scissors-operator shift. 

A perturbative electromagnetic field of frequency uj, and a wave vector q+G on a crystal 
produces a response of frequency uj and a wave vector q+G' (G and G' being reciprocal 
lattice vectors). The microscopic field of wave vector q+G' is produced by the umklapp 
processes as a result of the applied field Eq(q + G,uj) 

E (q +G,u) = J2 e G,G' (q, w)E(q + G\ uj) (4) 

G' 

where E(q + G, uj) is the total field which produces the non-diagonal elements in the micro- 
scopic dielectric function e^G^qj w). In the random phase approximation the microscopic 
dielectric function is given by |29| : 



e G ,G'(q, uj) = 6 G>G , - £ MX e^^-Lhs W' k + qi^ q+G)r k k) 

x (n, k|e- i ( q+G '>|n / , k + q) (5) 

Here n and n! are the band indexes, / n .k is the zero temperature Fermi distribution, and 
Q is the cell volume. The energies E n ^ and the the crystal wave function \n, k) are produced 
for each band index n and for each wave vector k in the Brillouin zone. 

The macroscopic dielectric function in the infinite wave length limit is given by the 
inversion of the microscopic dielectric function: 

e(uj) = lim — 1 — — — - — (6) 
[£G,G'(q» w )]o,o 

= e ,o(w) - lim q ^ Eg.gvo e o, G (q, uj)T^ G ,(q, ^)e G ',o(q, w) 

Where Tqq, is the inverse matrix of Tg,g' containing the elements cg,g' with G and 
G' 7^ 0. The first term of this equation is the interband contribution to the macroscopic 
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dielectric function and the second term represent the local-field correction to e. The most 
recent ab-initio pseudopotentials calculation found that the local-field effect reduces the 
static dielectric function by at most 5% |J. Previous calculations with the same method 
have also found a decrease of eoo by about the same percentage H,|T7| . We are looking at the 



effect of the local field using our all-electron basis-set; it should be of interest to compare 
all electron results with these obtained using the pseudopotential method. 

For insulators the dipole approximation of the imaginary part of the first term of equation 
(0) is given by @: 

e 2 (u) = ^-3- $3 / dk\(n,k\v\n' ,k)\ 2 f n;k (l - f n r ;k )d(e k)n , <n - fku) , (7) 

6uJ 71 n,n> J 

Here v is the velocity operator, and in the LDA v = p/m (p being the momentum operator), 
and where ek, n ,n' = E n i^ — E n>k . The matrix elements (nk|p|n'k) are calculated for each 
projection pj = ^dj, j = x or y and z, with the wave function |nk > expressed in terms 
of the full-potential LMTO crystal wave function described by equations ([I]) and @. The 
k-space integration is performed using the tetrahedron method |U| with 480 irreducible k 
points the whole Brillouin zone. The irreducible k-points are obtained from a shifted k-space 
grid from the high symmetry planes and T point by a half step in each of the k x , k y , and 
k z directions. This scheme produces highly accurate integration in the Brillouin zone by 
avoiding high symmetry points. 

To calculate these matrix elements we first defined a tensor operator of order one out of 
the momentum operator Vo = V 2 = and V±i = T^(-§^ ± The muffin-tin part 
of the momentum matrix elements is calculated using the commutator [V 2 , = 2V fl so 
that: 

J St 

St i*drM*£r + ^ +1) -/ y+1) )^(r) (8) 

where G x ^ nV m , are the usual Gaunt coefficients, and S T is the radius of the muffin-tin sphere 
of atom r. In the interstitial region the plane-wave representation of the wave function (see 
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equation ^) makes the calculation straightforward, but a special care has to be taken for 
the removal of the extra contribution in the muffin-tin spheres. However, we find it much 
easier and faster to transform the interstitial matrix elements as an integral over the surface 
of the muffin-tin spheres using the commutation relation of the momentum operator and the 
Hamiltonian in the interstitial region. The calculation of the interstitial momentum matrix 



elements is then similar to the calculation of the interstitial overlap matrix elements p3 
The k = case has been already derived by Chen using the Korringa, Kohn and Rostoker 
Green's-function method f32[ . We have tested that both the plane- wave summation and the 
surface integration provide the same results. 

Equation (|7|) can not be used directly to determine the optical properties of semiconduc- 
tors, when the GW approximation or the scissors operator is used to determine the electronic 
structure. The velocity operator should be obtained from the effective momentum operator 
p eff which is calculated using the self-energy operator, E(r, p), of the system |33 



v = p eff /m = p/m + <9£(r,p)/3p (9) 

GW calculations show that the quasiparticle wave function is almost equals to the LDA wave 
function [|T7],|TS|]. Based on this assumption, Del Sole and Girlanda show that the effective 
momentum operator p eff can be written in terms of the momentum operator p as follows 

(n, k|p eff \n, k) = (n, k|p|n, k)ejj£, tn /e k , n #, n , (10) 
where e^, n = E®, 1 ^ — E®^ is the difference between the quasiparticle energy E®, 1 ^ of the 



unoccupied state |n',k) and the occupied state |n, k). By substituting Equation |10| into 



equation [7], it can be easily shown |33j that in the case of the scissors operator, where all the 



empty states are shifted rigidly by a constant energy A, the imaginary part of the dielectric 
function is a simple energy shift of the LDA dielectric function towards the high energies 
by an amount A, i.e., e® p (uj) = ^^(u — A/h). The real part of the dielectric function is 
then obtained from the shifted €2 using Kramers-Kronig relations. The expression of ej*, is 
given by: 
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r QP_ i | 2e2 y n f \ |(n,k| P Kk)| 2 

6uj tt nn ,J (e kjn> + AJe ki<n 

e^ p is very similar to e^ )A except that one of the interband gap e ki „/ in is substituted by the 
QP interband gap e kjn / jn + A. 

To test for the accuracy of the calculation within the LDA the f-sum rule: 

2 EEWi-w^ = i, (i2) 



omn v k nn , ek,n',n 

where n v is the number of valence bands, is checked in all the calculations, and it is satisfied 
to within a few percent. 

It is easily seen that the dielectric function 62 calculated using the scissors-operator 
shift does not satisfy the sum rule (wp is the free-electron plasmon frequency): 

poo jy 

/ uje 2 {uj)duj = —Up (13) 
Jo 2 

because (i) e\ DA satisfies this rule, and (ii) e^ P is obtained by a simple shift of e\ DA by 
the scissors-operator A towards higher energies. Using the expression of the quasiparticle 
dielectric function in the scissors-operator shift approximation we show that e^ P satisfy the 
following integral sum rule: 

r^ P (u)du=^u' p 2 (14) 
jo 2 

where u' P 2 =u 2 P + £„,„' / dk\ (n, k|p|n', k) 1 2 /e kn> /„ ik (l - / <k ). We recover the usual 
sum rule when A is equal to zero. The non simultaneous satisfaction of both the f-sum 
rule and the integral sum rule given by Eq. 13 within the scissors approximation shows 
the limitation of this approximation. While the scissors operator approximation describes 
nicely the low lying excited states, which is seen in the good determination of the static 
dielectric function and the low energy structures, i.e. E 1 and E 2 , in the imaginary part 
of the dielectric function, it seems to fail for the description of the higher excited states. 
This is not surprising because the higher excited states which are free electrons like are most 
probably well described by LDA and need no scissors operator shift. This is supported by the 
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fact that the the energy-loss function, -Ime -1 , within the LDA has it maximum roughly at 
the free electron plasmon frequency whereas within the scissors approximation its maximum 
is shifted to higher energies as given by equation fll^]) . Fig. [I] show the energy- loss function 
of GaAs calculated within LDA (full curve) and within the scissors approximation (dashed 
curve). It is clearly seen that the maximum of the LDA curve has a maximum which is 
closer to the free valence electron plasma frequency of 15.5 eV. It is of general interest to see 
whether the calculated dielectric function within the GW approximation satisfies the integral 
sum rule. For our purpose the scissors-operator shift remains a good approximation for the 
description of the low-lying excited states of semiconductors and their optical properties. 



III. ELECTRONIC STRUCTURE OF SI, GE, AND GAAS 

The electronic structure of Si, Ge, and GaAs are obtained by solving the LDA equations 
by means of a full-potential LMTO basis-set as described above. Table I shows the orbitals 
used to describe the valence and conduction bands during the selfconsistency. The large 
number of orbitals used is found necessary to obtain converged excited states up to 5 Ry 
above the top of the valence states. However the total energy is insensitive to these high 
energy orbitals, but the presence of the 3d-core states of Ge and GaAs are important [ ^5] . 



Table II compares our band structure of Si for some high symmetry points with some re- 
cent results from first-principles calculations based on pseudopotential and Gaussian orbitals 



methods f27| , |28|| . We found a good agreement between our results and these calculations. 
This reflects the high accuracy of our unoccupied states which are used to determine the 
dynamical dielectric function. 

Table III shows the calculated equilibrium structural parameters, i.e., the electronic 
pressure and the bulk modulus at the experimental unit cell volume, Vq, and calculated 
cell volume, V. The calculated equilibrium volume V is at the most 2% smaller than the 
experimental value which correspond to a less than 0.5% deviation from the experimental 
lattice parameter. However the bulk modulus which is very sensitive to the slop of the 
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total energy versus the unit cell volume deviates at the most by 10% in the case of Ge, 
and when calculated at the experimental unit cell volume. But only by 5% when calculated 
at the theoretical equilibrium volume. Our calculation of the bulk modulus is in excellent 



agreement with other calculation pi ]. 

Figure ^ shows the LDA underestimated (a) direct band and (b) minimal gaps of Si, 
Ge, and GaAs compared with the Ge, and GaAs experimental results of Goni et al. |lj. For 
GaAs a cross over from direct band gap to indirect band gap takes palace between T and 
X at approximately 8 GPa. For Ge this cross over occurs along TL at a lower pressure of 3 
GPa. The direct band gap increases linearly with pressure and is in good agreement with 
the experimental results for both GaAs and Ge. There is no experimental data for Si under 
hydrostatic pressure. Table IV present the first and second order coefficients describing the 
dependence of the direct band gap at Y under hydrostatic pressure, E (P) = E + aP + bP 2 , 
compared to the experimental results of Goni et al. |l| . Apart for the underestimation of the 
band gap, the first and second coefficients of the pressure dependence of the band gap are in 
good agreement with the experimental results. This suggest that the scissors-operator shift 
is a good approximation for the description of the band gap under hydrostatic pressure. 



IV. OPTICAL PROPERTIES OF SI, GE, AND GAAS 

A. frequency dependent of the complex dielectric function of Si, Ge, and GaAs 

Figure and § present the imaginary part of the macroscopic dielectric function of 
Si, Ge, and GaAs obtained at the experimental ground state lattice parameters except for 
Ge where we have compressed the lattice parameter by about 1%. The compression is 
done because within LDA and at the experimental lattice parameter Ge is a semi-metal. 
The LDA e<i is shifted towards higher energy by the scissors-operator shift in order that the 
optical band gap agrees with experiment. The comparison to experimental results of Aspnes 
and Studna shows that all the features in the experimental spectra are reproduced by 
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the calculation. It is interesting to notice that the calculated LDA 62(00) of Si exhibit the 
largest underestimation of the E\ peak (about 50% in intensity) whereas, in Ge and GaAs 
the underestimation of the E\ peak is only about 12%. The E% peak is overestimated by 
LDA by about 34% for Si, 50% for Ge, and 60% for GaAs. This overestimation of the E2 
peak by LDA is due to a strong van Hove singularity near the X points of the Brillouin 
zone where parallel bands occur over a large plateau f5|,[^|. This overestimation can be 
reduced substantially by including the life-time broadening of the quasiparticles through a 
self-energy calculation. 

The effect of interband electronic transitions due to the 3<i semi-core states, without scis- 
sors operator shift, is presented in Fig. |5|a and refeps2highb for Ge, and GaAs respectively. 
For Ge the onset of transitions begun at photon energy of about 25 eV, and the intensity is 
very similar to the p-density of states of the empty states of Ge. This is because the 3d states 
of Ge are very narrow, and the dipole selection rules allow transitions only to the empty 
p states of Ge, the f-states in this energy range are absent. Whereas for GaAs, the onset 
of transitions begun at 12 eV, and the intensity €2 spectrum above 12 eV is very different 
from the empty p states of Ga. This is because of the relatively large dispersion of the 3c? 
semi-core states of Ga. It should be of interest to confirm experimentally these theoretical 
predictions. 

The real part 61(00) of the dielectric function of Si, Ge, and GaAs calculated by Kramers- 
Kronig transform of the imaginary part 62(00) are presented in Fig. ^ together with the 
experimental results of Aspnes and Studna || . In the same figure we have also presented the 
scissors-operator shift ef P (u;) and the high frequency asymptotic limit 61(00) = 1 — up/00 2 , 
where 00 p is the free-electron plasmon frequency. We notice that the analytic asymptotic 
limit matches nicely the calculated LDA 61, which is an indication of the quality of the cal- 
culation. For the P we need to use a different plasmon frequency as described in equation 
( due to the poor description of the higher excited states by the scissors approximation. 

In conclusion, we believe that the excitonic effects may be important for the dielectric 
function of Si but less for those of Ge, and GaAs. A QP calculation of the dielectric function 
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including the dynamical screening of the Coulomb interaction, like in the GW approximation 
of Hedin ||16|| , would certainly improve the intensity of at least the E 2 peak by introducing 
a life-time broadening of the quasi-particles. 



B. Hydrostatic pressure dependent of the static dielectric function of Si, Ge, and 

GaAs 

Figure [7] and || presents the hydrostatic pressure dependence of the static dielectric 
function, e^, of Si, Ge, and GaAs calculated within the LDA without and with the scissors- 
operator shift (SOS), respectively. Our data are compared to the experimental results of 
Goni et al. and to the pseudopotential calculations of Si and Ge of Levine and Allan || . 
Our calculation and the pseudopotential theory of Ref. || suggest that LDA is overestimating 
the static dielectric function of Si, Ge, and GaAs over the whole range of hydrostatic pressure, 
and that the use of a unique value of the scissors-operator shift for the correction of the band 
gap at T produces a nice agreement with the experimental results [|TJ. The static dielectric 
function decreases almost linearly with the pressure due to the increase of the direct band 
gap. However, for Si is almost constant with the pressure and this is because the increase 
of the direct band gap is almost compensated by a decrease of the indirect band gap (see 
Figure ||) . 

Table V, VI, and VII present the calculated pressure band gaps, static dielectric function 
and f-sum rule for Si, Ge, and GaAs, with a comparison to the experimental results of Ref. 

The agreement with the experimental results is excellent when the scissors-operator 
shift is used. The f-sum rule deviates at most by 5.2% from unity in the case of Ge which 
reflect the high precision of the calculation of the optical matrix elements. The fact that 
the f-sum rule is not quite exhausted for Ge and GaAs (deviation of about 5%) as compared 
to Si (deviation of about 1%) is not due to a possible incompleteness of our basis set |36| 



but rather to our use of all electron electronic structure. When the valence states are very 
well isolated from the core states, like in the case of Si where the core states lie about 
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80 eV below the valence bands, the sum rule should be exhausted. However, for Ge and 
GaAs where the semi-core 3c?-states are very close to the valence states and greatly affect 
the optical properties, the f-sum rule could deviate markedly from unity, i.e., the average 
effective number of electrons per atom contributing to the optical transitions is much larger 
than 4 electrons per atom |3"7|] . In pseudopotential theory, since the core states are absent, 
the f-sum rule is exhausted for all semiconductors ||. The details of the contribution of 
the 3d semi-core states to the oscillator strength and the study of the effective number of 
electrons contributing to the optical transitions are out of the scope of this paper and will 
be addressed elsewhere. 

The first and second-order coefficients describing the pressure dependence of the static 
dielectric function are presented in Table VIII. The results are compared to the experi- 
mental results of Goni et al. Ref. [l[], and the pseudopotential calculation of Ref. ||. The 
overall agreement with experiment and the pseudopotential calculation is excellent. 

In Table IX we present our calculation for the static dielectric function of Si, Ge, and 
GaAs including the spin-orbit coupling effect at the variational level, and the effect of the 3d 
states in the interband transitions. The calculated potential includes always the 3d states, 
and only the dielectric function is calculated with or without the 3d interband transitions. 
We have obtained that the inclusion of the 3d interband transitions increases slightly the 
static dielectric function, whereas the spin-orbit coupling increases it by 2.1% and 3.2% for 
Ge and GaAs, respectively. The of Si is insensitive to the spin-orbit coupling. The 
calculated scissors-operator shift 6^ including the spin-orbit coupling effect decreases by 
about 3.3% and 4.1% for Ge and GaAs, respectively. This is because the band gaps of Ge 
and GaAs are further reduced in presence of spin-orbit coupling which resulted in a larger 
scissors-operator shift for the determination of e^ p . 
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V. CONCLUSION 



The macroscopic dielectric function in the random-phase-approximation without local 
field effect has been implemented using the local density approximation with an all electron, 
full-potential linear muffin-tin orbital orbital basis-set. The method is used to calculate the 
optical properties of the semiconductors, Si, Ge, and GaAs, under hydrostatic pressure. We 
have found that the LDA overestimation the static dielectric function over all the pressure 
range from to 12 GPa, and that a single value of the so called scissors-operator shift which 
account for the correct band gap at T produced a good agreement with the experimental 
data of Goni and coworkers [[j]. This makes us conclude that because LDA underestimates 
the band gap it is incapable of producing the correct static dielectric function even though 

is a ground state property. 

Since (i) the Kohn-Sham (KS) density functional (DF) without the local density approx- 
imation should in principle produce the correct e^, and (ii) since the LDA calculation with 
the scissors-operator shift also produces the correct e^, we are tempted to conclude that the 
KS-DF theory should produce the correct band gap for semiconductors. This conclusion is 
not confirmed by a non-selfconsistent GW calculations which suggest that the true KS-DF 
theory also underestimates the band gap 

Our analysis of the dielectric function, the sum rules and the energy-loss function shows 
that while the scissors-operator shift is a good approximation for the low lying excited states 
it appears as bad approximation for the high energy excited states. This is because the high 
energy states are free electron like hence well described within LDA. 

Our calculation of the dynamical dielectric function shows that the E\ peak intensity 
is underestimated for Si by about 50%, and for Ge and GaAs by only 12%. These results 
imply that the excitonic effects may be important for the dielectric function of Si, but less 
for those for Ge, and GaAs. 

We have also shown that including the 3d semi-core states in the interband transitions 
hardly changes the static dielectric function, e OT , however their contribution to the intensity 
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of dynamical dielectric function for higher photon energies is substantial, and could be 
checked experimentally. We have also found that the spin-orbit coupling has a significant 
effect on eoo of Ge and GaAs, but not of Si. 

We thank J. W. Wilkins for interesting discussions. This research was supported in part 
by the U.S. Department of Energy Basic Energy Sciences, Division of Materials Sciences 
and by NSF, grant number DMR-9520319. Supercomputer time was provided by the Ohio 
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FIGURES 

FIG. 1. Calculated energy- loss function of GaAs within the LDA (full curve) and within the 
scissors approximation (dashed curve). It is clearly seen that the maximum of the LDA curve is 
much closer to the free valence electron plasma frequency of 15.5 eV. 

FIG. 2. Calculated (a) direct band gap Eo, and (b) minimum band gap E gap of Si, Ge, and 
GaAs as a function of hydrostatic pressure compared to the experimental results of Gohi et al. j!J 
for Ge (dashed line), and GaAs (thick line). Plot (a) shows that the direct band gaps increase 
almost linearly with pressure. Plot (b) shows that for GaAs there is a cross over of the band gap 
from direct to indirect at around 8 GPa, and a cross over for Ge at almost 3 GPa. The indirect 
band gap of Si decreases linearly with pressure. 

FIG. 3. Calculated imaginary part of the dielectric function of Si at the experimental equilib- 
rium volume, shifted by A = 0.6 eV towards higher photon energies, compared with the experi- 
mental results of Ref. j|. The experimental E\ structure at 4 eV is underestimated whereas the 
main E2 structure at 4.5 eV is overestimated. 

FIG. 4. Imaginary part of the dielectric function of Ge at 10 kbar shift by 0.4 eV and GaAs at 
the experimental equilibrium volume shifted by 1.1 eV, compared with the experimental results of 
Ref. |J. In both Ge and GaAs €2(0;) the experimental E\ is only slightly underestimated and E% 
is overestimated. 

FIG. 5. Contribution of the 3d interband transitions to the imaginary part of the dynamical 
dielectric function of (a) Ge (at 10 kbar) and (b) GaAs at the experimental equilibrium volume. 
The full line and the dashed line are with and without 3d interband transitions, respectively. Due 
to the narrow nature of the 3d semi-core states of Ge, the intensity of €2 above 25 eV is very 
similar to the empty p-density of states of Ge. Whereas for GaAs, the 3d semi-core states of Ga 
are relatively delocalized, which makes the intensity of €2 above 12 eV much different from the Ga 
empty p-density of states. 
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FIG. 6. Real part of the dielectric function of Ge (at 10 kbar) and GaAs, at the experimental 
equilibrium volume, compared with the experimental results of Ref. Q. The analytic asymptotic 
limit, shown by the empty circles, matches nicely the calculated spectra above 10 eV. 

FIG. 7. LDA scalar-relativistic calculated static dielectric function of Si, Ge, and GaAs as a 
function of hydrostatic pressure compared to the experimental results of Gohi et al. Q and the 
pseudopotential calculation of Levine and Allan ||] . 

FIG. 8. LDA plus the scissors-operator shift (SOS) calculated static dielectric function of Si, 
Ge, and GaAs as a function of hydrostatic pressure compared to the experimental results of Gohi 
et al. §. 
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TABLES 

TABLE I. Basis-sets used for the calculation of the excited states of Si, Ge and GaAs. Each 
orbital has different kinetic energy k 2 in its the interstitial region. For example, the 3s orbital of 
Si is used three times, each of the 3s wave function has a different kinetic energy in the interstitial 
region. 

Si: 3 x (3s, 3p) 2 x (3d) 

3x(4s,4p) 2x(4d,4/) 
Ge: 2 x (3d) 

3x(4s,4p) 2 x (id) 
2 x (5s, bp) 

GaAs: 2 x (Ga 3d) 

3x(Ga4s,4p) 2 x (Ga Ad) 

3 x (Ga 5s, bp) 
3x(As4s,4p) 2 x (As Ad) 
3 x (As 5s, bp) 
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TABLE II. Ei£ 


;envalues of Si at high symmetries 


points (r, X, 


and L) as compared to the 


results produced by 


means of a linear combination of gaussian orbitals 


p8|, and by pseudopotential 


(PP) method |7). ' 


1 he zero of energy is chosen at 1 25v 


i point. 




High-symmetry 


Gaussian 


PP 


Present 


point 


orbitals 




Calculation 




-11.91 


-11.91 


-11.96 


^25?;' 


0.0 


0.0 


n 

u 


Pise 


2.57 


2.55 


2.56 




3.24 


3.28 


3.20 


Xl v 


-7.77 


-7.76 


-7.82 


Xa 


-2.78 


-2.86 


—9 8^ 


X\ c 


0.65 


0.66 


0.62 




10.03 




10.03 


L-2v' 


-9.58 


-9.56 


-9.63 


Liv 


-6.94 


-6.96 


-6.99 




-1.17 


-1.20 


-1.19 


Lie 


1.47 


1.50 


1.44 


Lzc 


3.32 


3.33 


3.31 


L<2c' 


7.77 




7.66 


Indirect band gap 


0.52 




0.50 



25 



TABLE III. Calculated equilibrium volume (V), electronic pressure, and bulk modulus of Si, 
Ge, and GaAs. The the bulk modulus are calculated both at the experimental (Vo) and theoretical 
(V) unit cell volumes. The experimental results are shown between parenthesis. 



oemicoiiQucioi v o 


v / v o 




H{\ ) 




( A ) 






(LxFaJ 


(GFaJ 


Si 39 Q8 


QQO 


— 70 


95 8 












(98.8) 


Ge 45.27 


.988 


-.80 


71.0 


67.1 










(74.4) 


GaAs 45.12 


.984 


-1.2 


74.2 


69.3 










(74.7) 


TABLE IV. First and second-order coefficients 


describing the dependence of the direct band 


gap at T (Eq) under hydrostatic pressure (Eq(P) = 


= E + aP + bP 2 ) 


for Si, Ge, and GaAs. The 


experimental results are from Gohi et al. 


Ref. 0. 








Eo 


a (meV/GPa) 


b (meV/GP; 


a 2 ) 


Theory Expt. 


Theory 


Expt. 


Theory 


Expt. 


Si 3.273 


100.8 




0.05 




Ge -0.084 0.795 


125.4 


121 


0.2 


0.2 


GaAs 0.41 1.43 


99.1 


108 


-0.1 


-0.1 
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TABLE V. Calculated pressure, band gap, static dielectric function with and without scis- 



sors-operator shift (SOS), and the f-sum rule of Si as 


a function 


of volume. pq|, 




v„ /V 
Vo/ v 






T Oil Y~V~1 


gap^ev ) 


£oo (LDA) 


eoo (SOS) 


eoo (Expt.) 


1.000 




-.75 


0.988 


0.50 


13.75 


12.08 


12.0 a 


1.025 




1.8 


0.989 


0.46 


13.65 


12.00 




1.050 




4.8 


0.987 


0.42 


13.61 


11.96 




1.100 




9.8 


0.991 


0.34 


12.57 


11.98 




a Ref. 


38 














TABLE VI. Calculated pressure 


, band gap, static dielectric function with and without scis- 


sors-operator shift (SOS), and the f- 


sum rule of Ge 


as a function of volume. The experimental 


data are from Gohi et al. 


0- 










V /V 




P(GPa) 


f-sum 


gap(eV) 


eoo (LDA) 


eoo (SOS) 


eoo (Expt.) 


1.025 




1.0 


1.053 


.04 


19.24 


15.32 


15.59 


1.050 




2.9 


1.041 


.21 


18.14 


14.71 


15.01 


1.075 




5.0 


1.042 


.29 


17.44 


14.33 


14.49 


1.100 




7.2 


1.043 


.37 


16.70 


13.90 


14.07 


1.150 




12.3 


1.045 


.48 


15.81 


13.46 


13.63 
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TABLE VII. Calculated pressure, band gap, static dielectric function with and without scis- 
sors-operator shift (SOS), and the f-sum rule of GaAs as a function of volume. The experimental 
data are from Gohi et al. Q. 



V /V 


P (GPa) 


f-sum ; 


gap(eV) 


£oo (LDA) 


£oo (SOS) 


(Expt.) 


1.000 


-1.2 




1.041 


.29 


14.44 


11.0 


11.05 


1.025 


.68 




1.042 


.48 


13.93 


10.72 


10.88 


1.050 


2.8 




1.043 


.66 


13.45 


10.53 


10.69 


1.075 


4.7 




1.044 


.85 


13.09 


10.41 


10.53 


1.100 


6.8 




1.044 


1.02 


12.75 


10.25 


10.34 


1.150 


11.8 




1.046 


1.08 


12.20 


10.03 


9.90 


TABLE VIII. 


First and second-order coefficients describing the dependence 


of the static di- 


electric function on hydrostatic pressure (e 0O (P) 


= e oo + aP + bP 2 ) 


for Si, Ge, 


and GaAs. The 


experimental data 


are from Gohi et al. 












e° 


i 


a(l/GPa) 


b (1/GPa 2 ) 


dln( eoo )/dP(10- 12 /Pa) 




Theory 


Expt. 


Theory 


Expt. 


Theory Expt. 


Theory 


Expt. 


Si 


12.05 




-0.032 




0.0025 


-2.65 




Si 


11.16* 




-0.027 a 




0.0013 a 


-2.6,-2.43 a 




Ge 


15.58 


15.94 


-0.32 


-0.36 


0.012 0.014 


-20.21 


-22.60 




16.04 a 




-0.46 a 




0.018 a 


-31,28.66 a 




GaAs 


10.83 


10.92 


-0.11 


-0.09 


0.004 


-10.43 


-8.06 



a Pseudopotential calculation of Ref. ||, slightly larger numbers are quoted for dln^oo) / dP 
in their Table VI. 
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TABLE IX. Calculated static dielectric function of Si, Ge, and GaAs at the equilibrium lattice 
parameter (except for Ge where it is calculated at a slightly smaller lattice parameter (1% smaller) 
than the experimental one because Ge is a metal in LDA for V/Vq = 1.). The calculation are 
done using scalar relativistic (SR) LMTO without 3d states, with the 3d states (SR+3c£), with the 
spin-orbit coupling at the variational level (SR+SO), and with the SO coupling and the 3d states 
included (SR+SO+3d). 







Si 




Ge 




GaAs 




LDA 


LDA+SOS. 


LDA 


LDA+SOS 


LDA 


LDA+SOS 


SR 


13.75 


12.08 


18.14 


14.71 


14.44 


11.0 


SR+3d 






18.16 


14.73 


14.47 


11.03 


SR+SO 


13.69 


12.0 


18.52 


14.23 


14.90 


10.52 


SR+SO+3<i 






18.54 


14.25 


14.93 


10.55 


Expt. 




12. a 11.4 b 




14.98 c 




10. 9 a 



Ref. BHl 



Ref. 133 



Ref. § 
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